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We use the (1 + l)-dimensional Kardar-Parisi-Zhang equation driven by a Gaussian white noise 
and employ the dynamic renormalization-group of Yakhot and Orszag without rescaling [J. Sci. 
Comput. 1 , 3 (1986)]. Hence we calculate the second and third order moments of height distribution 
using the diagrammatic method in the large scale and long time limits. The moments so calculated 
lead to the value S = 0.3237 for the skewness. This value is comparable with numerical and 
experimental estimates. 
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I. INTRODUCTION 

The study of surface growth has been one of the most important problems in non-equilibrium statistical physics 
over the past few decades [IHl- The most generic continuum model of surface growth is the Kardar-Parisi-Zhang 
^PZ) equation that is endowed with interesting properties of statistical scale invariance. Kardar, Parisi and Zhang 
suggested a nonlinear differential equation for local surface growth in the form 

y(V/i)2-|-? 7 (x,t), (1) 

where /i(x, t) is the height of the surface at position x and time t on a d dimensional substrate, vq is the surface 
tension that has a tendency to make the surface smooth, and the coupling constant Ao measures the strength of the 
nonlinear interaction term. The nonlinear term induces local growth along the normal to the surface and gives rise 
to lateral correlations. On the other hand, the linear term (containing z/g) is responsible for diffusion of particles to 
the local minima Q. The driving term r]{x,t), describing the random deposition of particles, is assumed to obey a 
Gaussian distribution to account for the stochastic nature of the flux of particles. It is taken to be a Gaussian white 
noise with zero mean, t)) = 0, and with correlation 

(? 7 (x, t) r]{x', t')) = 2Do - x') 6{t - t'), (2) 

where Dq is a constant and the angular brackets denote ensemble averages. 

There are many deposition models that have been identified with the KPZ universality class. A few examples are 
the ballistic deposition d, Q , the Eden model [iSIl , the restricted solid-on-solid (RSOS) model 0, and the single step 
model (SSM) [1, [l3|- A large number of growth experiments show scaling exponents close to those of the KPZ growth 
problem. A few important phenomena are thin film growth [I| , bacterial colony growth [IS El, growth of fractals [11 , 
turbulent liquid crystal (TLC) growth [IS El, and one dimensional polynuclear growth (PNG) [TtI - I^ . etc. Apart 
from such growth models, the KPZ problem is related to various other processes such as the noisy Burgers equation 
[^ , flame front propagation (^ . , directed polymer in random media [liHiil, interface roughening due to 

impurities (ttI . , and growing interfaces in randomly stirred fluids [^ . A great amount of work has been carried 
out, mostly via numerical and experimental studies, in various KPZ type surface growth problems. The interplay of 
non-linearity, surface tension, and uncorrelated noise in such problems establish a universality class distinct from that 
of the Edward-Wilkinson type growth in the large-scale long-time limit. The mean-square of the height fluctuations 
is related to the critical exponents [s^ as 

([/i(x, t) - h{x\ - |x - x'px ^ ^ ( 3 ) 
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where x is the roughness exponent describing the self-afHne geometry of the surface, z is the dynamic exponent (the 
ratio j = Bis the growth exponent), and is a scaling function. The roughness exponent x is important 
parameter in the studies of adsorption, catalysis (3^ , and optical properties [3^ of a thin film. The properties of 
a rough surface are determined by the distribution of height fluctuations and it deserves attention both in theoretical 
and experimental studies of growing interfaces [^. 

Various analytical approaches have been employed to study the universality class of the KPZ equation on the basis 
of scaling exponents in different dimensions. The dynamic renormalization-group by Kardar, Parisi, and Zhang Q 
leads to the values of roughness exponent X — ^ and dynamic exponent z = | at one-loop order for the (1-1-1) 
dimensional KPZ equation. Motivations in the theoretical study of the KPZ equation in higher dimensions have led 
to formulations of different analytical techniques. Examples of such theoretical studies are the mode coupling scheme 
( 334 ^ . th e operator product expansion , the self-consistent expansion [s^ and a nonperturbative renormalization 
group [39l . |40| for the calculation of scaling exponents in the strong coupling regime. 

These exponents have also been computed numerically considering different growth rules. Apart from the numerical 
studies, many experiments have been carried out to find these exponents. Various experimental studies BEllil 
have indicated that the roughness exponent is about 0.50 and the growth exponent 0.33 which have been identified 
with the universality class of the (1 -I- 1) dimensional KPZ type surface growth. 

Besides the critical exponents, the probability distribution function is an important feature to classify the universal¬ 
ity of a physical process [ 2 ^. In experiments, measurements of normalized moments is expected to be more accurate 
than the measurement of scaling exponen ts |4| . Thus higher order moments can infer about the universality class in 
a better way than the critical exponents [43l| . Consequently, higher order moments are very important in the study 
of surface morphology and the universality class can be better realized through the values of higher order moments 
and a parameter such as skewness. 

The skewness has been computed numerically employing a variety of deposition algorithms. Krug et al. (43 |. using 
the simulation in a single step model for flat initial condition, obtained [S'! = 0.28 ± 0.04 in the transient regime. 
Following the same model, they prep ared stationary interfaces by taking uncorrelated spins (ct^ = ±1) and obtained 
~ 0.33. Prahofer and Spohn [20j took the polynuclear growth model and mapped it into a random permutation 
through the droplet geometry thereby onto Gaussian random matrices to understand the dependence of the initial 
conditions on height fluctuations. They inferred that the droplet and flat substrates have the same scaling form but 
distinct universal distributions. They estimated the skewness for three different shapes, namely, curved, flat, and 
stationary self-similar in (1 + 1) dimensions. For the flat shape, they obtained S = 0.2935, for the curved shape, 
S = 0.2241, and for the stationary self-similar case, S = 0.35941. They proposed an expression for the height 
distribution, namely h{x,t) ~ Voot + (rt)^/^C with C a random variable, where P = is a model parameter and 

Voo is the growth rate in the asymptotic limit [l^. It was found that C obeys the Tracy-Widom (TW) distribution 
corresponding to the largest eigenvalues of random matrices [l^ . For curved interfaces the random matrices form a 
Gaussian unitary ensemble (GUE) whereas for flat interfaces they form a Gaussian orthogonal ensemble (GOE). 

In an experiment on growing interfaces in liquid crystal turbulence, Takeuchi et al. [3 [3 found that the growth 
and roughness exponents are the same as those of the KPZ type growth in one dimension in the asymptotic limit. 
Their experimental data indicated the value for skewness S ~ 0.29 for a flat interface, whereas for a curved interface 
their experimental data converged to ^ = 0.2241. They concluded that the probability distribution function (pdf) 
of interface fluctuations precisely agrees with the GOE of TW distribution for the flat interface, whereas the curved 
interface fluctuations agree with the GUE of TW distribution, up to fourth order cumulants. Sasamoto and Spohn 
[313 solved the (1 -I- 1) dimensional KPZ problem with an initial condition of curved-height profile and showed that 
the pdf follows the GUE of TW distribution of random matrices. 

It may be noted that there have been very few analytical evaluations of the skewness and higher order moments 
for the KPZ type growth problem. The one known to the authors is a mean field calculation yielding S = ±0.46 in 
(1 ± 1) dimensions [3 with the flat initial condition h{x, 0) = 0 for the transient regime. 

In this work, we are interested in the KPZ growth problem for a flat interface and seek to calculate the skewness 
of height fluctuations in the stationary state. Gonsequently, we apply the dynamic renormalization group scheme 
without rescaling to the KPZ equation. This scheme was previously employed by Yakhot and Orszag |48l | to calculate 
various universal numbers in the case of hydrodynamic turbulence. This scheme enables us to calculate the second and 
third order moments of height fluctuations in a straightforward manner. The ensuing result for skewness is compared 
with the findings of various numerical, experimental, and theoretical studies in Table 1. 

The paper is organized as follows. In Section 2, the renormalizaiton-group scheme without rescaling is applied 
to the KPZ problem. Section 3 outlines the definition of statistical moments of height fluctuations and presents 
calculations of the second and third order statistical moments. Finally, Section 4 presents a discussion and conclusion 
and a comparison with other findings. 
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II. RENORMALIZATION SCHEME WITHOUT RESCALING 

The nonlinear dynamics described by the KPZ equation ([1]) incorporates interaction among many degrees of freedom 
[33 • The complexity of such interactions among the collective set of height fluctuations is most easily seen when we 
Fourier transform the height fluctuations /i(x, t) and the driving field r](x,t). The Fourier space is also suitable for 
employing the dynamic renormalization-group techniques [i^ . The Fourier transform of the height fluctuations /i(x, t) 
is expressed as 

= (4) 

where d is the substrate dimension. The stochastic noise is also Fourier transformed in a similar manner. The 

Fourier amplitude of the noise fluctuations has a zero mean, (? 7 (k, ui)) = 0, and the noise-correlation can be expressed 
as 


(r 7 (k, uj) rjCk',io')) = 2 Dq (27r)'^i5‘^(k -f k') 2tt5{uj + w'), (5) 

in the Fourier space, as a consequence of Eq. [51 Using Eq. (jT]), the Fourier transform of the KPZ equation [Eq. (|T])] 
is obtained as 

(-■iw -f v^k^) /i(k, w) = ??(k, w) - Y JJ ^ 2 ^)^+! “ q)] - q, w - fl). (6) 

which is in a form particularly useful for implementing the renormalization-group scheme. 


A. Scale Elimination 

To implement the renormalization-group scheme, we eliminate height fluctuations (q, fl) belonging to the shell 
Ape”’’ < 9 < Aq in the wavevector space by substituting for /i^(q, fl) in the equation for /i^(k, w) following from 
Eq. (jb)). This process generates a perturbation series in powers of the coupling constant Aq. Considering terms up to 
second order in Aq yields the equation 

[—iuj + h'ok'^ + S(k, w)] (k, uj) = 77^ (k, uj) 

-y J j (k-q)]/i<(q,f^)/i<(k-q,a;-f2). ( 7 ) 

in the range 0 < fc < Aoe”’’ in the wavevector space, where Il(k, w) is the self energy correction represented by the 
amputated part of the Feynman diagram shown in Fig. 1. 



FIG. 1: Self-energy correction. The self-energy S(k, w) corresponds to the loop. Propagators are indicated by arrowed lines 
and correlation by a wiggly line. 


The corresponding loop integral is given by 

2 ^ jd. 


S(k,w) = 4 - 


Ao 


(27r)‘^ 


(k • q) [q • (k - q)] 


27r 


|G>(g)|2(2Uo)G>(fc-g), 


( 8 ) 


where Go{k) = Gnfk, ui) = [—ioj + VQk'^]~^ is the bare propagator and the prefactor 4 is a combinatorial factor. 
Following Refs. |30l.l48l|. we symmetrize the internal momenta by taking the transformation q —>■ (q-|-k/2). Performing 
the frequency convolution and evaluating the integral over the internal momenta in the shell Ape”’’ < g < Aq yields 
the self energy 


AgUp Sd (2-d\ 

2r'gAg“‘^ (27r)'^ 


1 


e 


E(/c,0) = 


2d 


2-d 


(9) 
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in the large scale {k —>• 0) and long time (w —)■ 0) limits, where Sd = r(d/ 2 ) surface area of a sphere of unit radius 

embedded in a d dimensional space. As a result of the above elimination, the effective surface tension is obtained as 


i/< (r) = Vo 


Ik 

4: 

VqI^-o 


,{2-d)r _ I 


( 10 ) 


where Kd = and the second term in the parentheses comes from the self energy correction. 

The height-height correlation is also expanded in a perturbative series in a similar manner. This gives rise to a 
correction to the noise amplitude, given by 


2D<{r) = 2Do + 2 



dS 

(27r)'^ 


[q.(k 


q)]^ S i^Do)^ |Go>(fc - <?)P, 


( 11 ) 


where D^{r) is the effective amplitude of noise correlation, whereas Dq is the bare parameter appearing in the noise 
correlation in Eq. ©■ The corresponding equation is shown diagrammatically in Fig. 2. 




FIG. 2: Perturbation expansion of the correlation Q(k,a;) to one-loop order. 


Calculating the loop integral in the large scale and long time limits, fc —?> 0 and w —>■ 0, the correction to the noise 
amplitude is obtained as 


AD = D 

2-d ■ 

Thus the effective amplitude of the noise correlation is given by 


D<{r) = Do 


1 \lDo 

^ 2-d 


( 12 ) 


(13) 


We observe that the surface tension vq and noise amplitude Dq acquire corrections due to the elimination of small 
scales belonging to the high-momentum shell Aoe~^ < fc < Aq. 


B. Flow Equations and Fixed Point 


To implement the renormalization scheme, we shall follow a procedure suggested by Yakhot and Orszag [3 [13| 
where the renormalized parameters are not rescaled after the above scale elimination operation. A particular advantage 
with this scheme is that the flow equations for the renormalized parameters are obtained directly with respect to the 
elimination parameter r. Implementing this scheme, we obtain, from Eqs. [TU] and 1131 the flow equations for the 
renormalized surface tension v{r) and renormalized noise amplitude D(r) as the differential equations 


and 


^ = / 2-d \ XlDjr) 

dr 4 \ d J v'^(r)A‘^~'^(r) 


1 \lD\r) 

dr 4 v^{r)A‘^~'^{r)’ 


(14) 


(15) 


where A(r) = Ape In this scheme, there is no flow equation for the coupling constant Aq as it does not acquire any 
correction due to Galilean invariance. In order to find the fixed point, we define an effective coupling, g(r), as 


g{r) = Kd 


AgJ(r-) 

v^{r)A^~'^{r) 


(16) 
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Using Eos. [HlandfTSl the flow equation for this effective coupling is obtained as 

^ = ag{r)-bg^{r), (17) 

where a = 2 — d, and b = Integrating this equation, we obtain an r-dependent expression for the effective 

coupling, given by 


gir) = 


50 


1 + !5o(e“’- - 1 ) 


(18) 


where go = 5 ( 0 ) = Kd 3 %-ti - The fixed point value g* is obtained in the limit r —>■ oo. For d < 2, we get 


’'o^o 


a 2d{2 — d) 
b ^ (3 - 2d) ■ 


(19) 


We see that the fixed point value g* diverges for the substrate dimension d = 1.5 and it is finite and positive in the 
range 0 < d < 1.5. However, in the range 1.5 < d < 2, the coupling constant is finite but ne gat ive, and it vanishes 
at d = 2. These fixed point values are consistent with Frey and Tauber’s one-loop calculation [^, Cf. Eq. (3.18)]. In 
this paper, we are interested in the substrate dimension d = 1 ; thus the effective coupling constant approaches the 
fixed point value g* = 2. 

Using Eqs. (USD and CHI), the differential equations d and C3 yield the exact solutions 


and 


u{r) = Vo 


l + ^(e“’--l) 
a 


aj ^hd 


D{t) = Do 


l + ^(e“’--l) 
a 


1/46 


( 20 ) 


( 21 ) 


For very large r, the above solutions lead to the asymptotic expressions 


and 


/ \ f ^50 ar 

v(r) ~ Vo -e 


D(r) ~ Do i 


aj Ahd 


1/46 


( 22 ) 


(23) 


in the large scale limit. Noting that a = 1 and 6=5 ^r our case d = 1, these expressions for surface tension and 
noise amplitude reduce to 


^ r/2 


( \ / ‘ 
v{r) ~ vo\l—e 


(24) 


and 




(25) 


These asymptotic expressions, for very large r, correspond to the renormalized surface tension 




■^ 0-^0 ^-1/2 


(26) 


and renormalized noise amplitude 


D{k) ~ Doi 


'AqUo ,,- 1/2 
27rz.3 


(27) 
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in the large scale long time limit. 

The dynamic exponent z can be defined via the renormalized response function as 

G~^{k,oj) = \—iuj + v{k)k'^'\ ^ oc k^ (j) , (28) 

suggesting the scaling v{k)k‘^ ^ k^. This leads to the dynamic exponent 2 = f and roughness exponent X = ^, the 
latter being a consequence of the scaling relation x + z = 2 . 

III. STATISTICAL MOMENTS AND SKEWNESS 

The nth moment of the height fluctuations is defined as 

Wr, = {[h{^,t)-h{t)r). (29) 

These moments obey power laws in the stationary state and they scale as Wn ^ where L is the size of the 
substrate. 

The statistical measure corresponding to the (square of) interface width (or standard deviation) is given by the 
second moment 


W 2 = - (h(x,t)}^. 

The skewness is related to the third moment 


(30) 


W 3 = (h^(x, t)) - 3(h^(x, t)){h{x, t)) + 2(/i(x, t))^. 


(31) 


In this paper, we calculate the skewness S of surface height fluctuations in the KPZ surface growth model. It is 
defined as 


W 3 

^ ^ (IT2)3/2 ■ 

We present the calculations of the moments W 2 and W 3 in the following subsections. 


(32) 


A. The Second Moment 


The second moment is expressed in the Fourier space as 


{h‘^ix,t)) = 


dfikdw fdfik'duj' 


{2ttY+^ J (27r)^+i 


(/i(k, w) h(k',a;')) e 


' p*(k+k')'X 


(33) 


We shall assume the growth process to be statistically homogeneous in space and stationary in time. This assumption 
yields the form 


(/i(k, uj) /i(k', uj')) = Q(k, oj) {2 ttY ^“(k + k') 27r (5(a; + w') 


(34) 


From Eq. we see that (/i(k, w)) = 0 for any k 7 ^ 0, implying {h{x,t)) — 0 for all practical purposes. Thus from 
Eqs. [Sni [Ml and [Ml we obtain 


X,. /. 9 / f (Ykdoj , 

W 2 = {h (x,t)) = J Q(k,a;) 


We write the integrand in terms of renormalized quantities as 

d'^kdbJ 


W 2 = 


J { 2 nY+^ G(k,a;)£ 2 (k,a;)G(-k,-a;) 


(35) 


(36) 


We hrst consider the bare value 


Lf(k,..)=2Do + 2 (^)'y' [q.(k-q)p|Go(q,fI)P |Go(k - q,.. - fI)p(2Do)^ (37) 
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where the propagators are unrenormalized. 

We evaluate the second term in Eq. [33 which corresponds to the amputated part of the loop diagram in Fig. 2. 
Performing the integrations over the internal frequency and internal momentum in the shell Age”’' < q < Ag, we 
obtain 


L<{r)=2Do + Kd 




2i/3Ag-^ 2 - d 


(38) 


Following Yakhot and Orszag’s procedure of renormalization, we make the assumption that thin shells in momentum 
space are eliminated recursively in iterative steps. This leads to a differential equation for L 2 (r), 


^ 1 XID\t) 
dr 2 tt {r) A{r) ’ 

representing the evolution of L 2 {r) with respect to the recursive steps of the shell elimination scheme. Using Fas. [Ml 
and [2^ and integrating over r, Eq. |39| yields 


L2{r) = Do 


\2\lDo ^,12 


ttz/Mi 




(40) 


for d = 1, in the asymptotic limit of large r. We transform this expression into a wavenumber and frequency dependent 
expression identifying Age”’’ as kf{-^) where z is the dynamic exponent and /(•) is a dimensionless scaling function. 
Thus, we obtain the renormalized function corresponding to Eg. 1401 as 


L2(k, Lo) — Do 




(41) 


We identify the scaling function by considering consistency in the w —>■ 0 limit, so that 


k 



1 

v'^{k) |G'(k, a;)p 


(42) 


where the modulus of the response function G'(k, ui) signifies further consistency with the fact that Age ’’ is a real 
quantity. Thus the renormalized quantity L 2 (k, w) is expressed as 

L2(k,cc) = ^fc|G(k,cc)|. (43) 

We notice that the first diagram in Fig. 2 does not contribute to L 2 and the contribution comes solely from the loop 
diagram. 

Substituting the expression |43| in Eq. |36| and treating the propagators as renormalized given by Eq. |28l with the 
renormalized surface tension I'ik) coming from Eq. 1261 we obtain the contribution to the second moment as 


\2 p )2 

W 2 = ^ 


Performing the frequency integration using 


f*+oo 


d^k r+“ ^ 

( 27r )'^ 2tt 


13/2 


+ iy‘^{k) k'^ 


du) 


_ r(a - \) 

J-00 (^^ + (m^)“”^/^ r(a) 

and carrying out the momentum integration in Eq. 1441 we obtain 

w 4: f Do \ I 
W 2 = 


TT \ 27ri/g J ^ 


(44) 


(45) 


(46) 


where /i is an infrared cutoff in the momentum integration. 


















B. The Third Moment and Skewness 


The third moment t)) can be expressed in the Fourier space as 



(/i(k, oj) h{k', uj') /i(k", uj")) e*(k+k'+k")-x 


Contribution to W 3 comes from the one-loop diagram shown in Fig. 3. 
one-loop contribution L^(k, k') as 


d’^k" duj” 
(27r)‘^+i 


(47) 


Thus W 3 can be written in terms of the 


n fld+lp ^ ^ ^ ^ ^ ^ 

( 2 ^ J ( 2 ^ 

where k stands for (k, w) and k' for (k',a;'). 

We first consider the bare value of the loop integral. Integrating over oj”, the bare loop integral can be written as 

4 °^(k,a;;k',a;') = 8 J [(q - k) • (k'-f k - q)] [q • (k - q)] 

[-q • (q - k' - k)] Qo(q, (5o(k - q,a; - S2) (5o(k + k' - q, w -I- w' - fl) (49) 



FIG. 3: The third-order moment. 


Carrying out the frequency convolution in Eq. 1491 we extract the leading order contribution from this integral in 
the large scale and long time limits, namely the limits A: —^ 0, A:' —>■ 0, w —^ 0, w' —>■ 0 for the external momenta and 
frequencies. Working out the momentum integration in the high-momentum shell Age”’' < 9 < Ag, we obtain 




3 n3 p(4-d)r _ 1 

2 4-d • 


(50) 


As before, we consider the iterative nature of the shell elimination scheme in thin shells in the momentum space, and 
obtain the flow of L^ir) in the form of a differential equation 


dLs _ 3 X^D^ir) 1 
dr 2 tt v^{r) A^{r) 


for d = 1. The functions z/(r) and D{r) being known from Eas.[24land[2^ the differential equation is solved to obtain 
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in the asymptotic limit of large r. 

The corresponding renormalized function L^(k]k'), being symmetric with respect to interchange of k and k', its 
frequency dependent expression can be obtained by replacing (Aoe“’’)“^ with the expression 




Employing Eq. 021 in 1121 we thus obtain 


Q \5 r)4 

L3(k,cc;k',a;') = ^^ |G(k,a;)p \G{\^',J)\\ 


Using this expression in Eq. 1481 the third moment is obtained as 
3 f dfi+^k f d^+^k' 


W, = 


Stt^ Vq J (27r)‘^+^ J (27r)'^+^ 


k^k'^ G{k) |G(fc)r G{k') |G(fc')r Gi-k - k') 


(53) 


(54) 


(55) 


where the propagators are treated as renormalized as expressed in Eg. I28l with the renormalized surface tension v{k) 
given by Eq. 1261 Performing the frequency integrations over oj and oj' yields 


0/7^ \ 3/2 /■ + 00 P + OQ 


in one dimension, where 


F{k^,k'^) = 


UiK,K) 

y{kx,K) 


with 


and 


V{k,X) = 16 |A:.| Ifc^l (|fc,|3/2 + |fc;|3/2 + \k, + fc;|3/2)' . 


The integrations in Eq. [56l can be decomposed to obtain 

3/2 


o / ^ \ o/^ • I>oo poo poo poo 

rJ "^^“7 dk',Fik,,k',)+2j dk^j dk',F{-k,X) ■ 

where we have set infrared cut offs at ^ as these integrals have infrared divergences. Thus we write 

o / n \ 3/2 

PCM) + 2j(rti 


where 


and 


poo poo 

I{k)= / dkx / dk'^F{kx,k'^) 

J U U 


pOO pOO 

/ dk^ / dk'^F{-kx,k'^). 

J pL J pi 

The infrared divergences in these integrals suggest the following forms 

/(m) = /oA^-'/" 

J(Ai) = Jo/i-"/" 


(56) 


(57) 


C/(fc,, fci) = 3 (|fc,|3 + |fc;|3) + 4|fc, + {\k4^G + |fc;|3/2^ + lA\kxf^\k'^^/^ + \kx + k'S (58) 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 

(65) 
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where /q and Jo are dimensionless constants. 

Substituting from Eqs. (m and (1651) in Eq. dSH), we obtain the third moment as 

/ D \ 1 

= 3 f ^ ) (/o + Jo) 


2TTUn 




3/2 


According to the definition of skewness, we thus obtain from Eqs. |4^ and [66] 

We calculate the constants lo and Jo from Eqs. |62]and|63l using the expressions given by Eqs. 
obtain 


( 66 ) 

(67) 

and mi and 


lim = 0.034946 

At-»-0+ 

(68) 

lim [/r3/2j(Ai)] = 0.120089 

(69) 


by means of numerical integrations. The computation shows convergence to the above numerical values as the 
parameter /r is chosen to approach values close to zero. 

Erom Eq. Ezl the value of skewness is thus found to be 

S' = ^ (0.155035) 7r3/2 = 0.323732. (70) 


IV. DISCUSSION AND CONCLUSION 

We employed Yakhot and Orszag’s scheme of renormalization without rescaling and obtained the renormalized 
surface tension and the strength of the noise correlation for the surface growth problem governed by the KPZ dynamics 
on a flat substrate. This scheme of renormalization is slightly different from the usual perturbative renormalization 
group analysis with rescaling that has been employed for dynamical problems by Ma and Mazenko [5^ , Eorster et al. 
[2ll] . and Medina et al. [^. This method allowed us the advantage of obtaining the flow equations directly without 
rescaling by considering the iterative nature of the scale elimination procedure. This yielded the fixed point from the 
r-dependent expression of effective coupling constant g{r) in the limit r —> oo. Similar to the other calculations, the 
renormalized surface tension and the strength of the noise correlation are found to be renormalized in the same way so 
that D(r)/vlr) is r-independent, a consequence of fluctuation dissipation theorem for the case of (1 + 1) dimensional 
KPZ equation 

To obtain a numerical value for the skewness, we employed the diagrammatic approach for the second and third 
order moments W 2 and W3. The Eourier integrals of these moments involve the loop integrals L 2 and L3. The 
simplicity of Yakhot and Orszag’s renormalization scheme allowed us to find renormalized expressions for these loop 
integrals in a straightforward manner. Although the renormalized diagrams are infrared divergent, the calculated 
value of skewness turns out to be finite due to cancellation of the infrared cutoff parameter /r. We obtained a value of 
skewness S = 0.323732 for the flat geometry in the stationary state which is compared with the results of numerical 
simulations for various growth models and those of experiments in Table 1. We present a discussion with regard to 
these results in the following paragraphs. 

It has been shown by numerical simulations for polynuclear growth (PNG) that the roughness and growth exponents 
are in good agreement with the one dimensional KPZ exponents [l^Il8j|. Eurther numerical work by Krug et al. [43 | 
and Bartelt and Evance [IJ have ensured that the PNG model belongs to the universality class of the KPZ growth 
model 0. Prahofer and Spohn [ 2 ^ have shown that the PNG model follows the TW distribution with different initial 
conditions. They estimated the skewness for three different shapes, namely, S = 0.2241 for the curved shape (GUE 
TW), S = 0.2935 for the flat shape (GOE TW), and S = 0.35941 for the stationary self-similar case [l^. On the 
other hand, the distribution of height fluctuations for the KPZ growth model with sharp wedge initial condition was 
shown to be the same as that of the GUE TW distribution, as established by Sasamoto and Spohn [4^, |4^ . Galabrese 
and Doussal obtained the GOE TW distribution by mapping the one dimensional KPZ problem with flat initial 
condition to a one end free directed polymer, referred to as a point-to-line configuration. The curved initial condition, 
on the other hand, maps on to a point-to-point configuration of the directed polymer. 

For the case of directed polymers at zero temperature in a random potential (DPRP), Kim et al. [ 2 ^ introduced two 
types of random site potentials /i(x, t), namely, uniform and Gaussian distributions for /x(x, t), with the bending energy 
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( 7 ) of the polymer as the only tunable parameter. They obtained skewness S = —0.29 ± 0.02 of the minimum energy 
distribution in 1 + 1 dimensions for uniform distribution of /i(x, t) for a point-to-line configuration via simulations for 
7 > 1. The same value of skewness was obtained for Gaussian distribution of /r(x, <), which is independent of 7 . Kim 
and others [ 2 ^, studied height fluctuations of surface growth using the RSOS model with a flat initial condition 
where the scaling form of the height distribution matches with the energy distribution of the DPRP within numerical 
accuracy. The skewness in the same model turned out to be 5" « —0.29, suggesting universality of the probability 
distribution function. 

Using a mean field theory in terms of densities at different heights applied to the KPZ equation in (1 + 1) dimensions, 
Ginelli and Hinrichsen started with the flat initial condition h{x, 0) = 0 and obtained the skewness S = ±0.46 for 
the transient regime. Takeuchi et al. [l^ carried out an experiment on a growing interface in liquid crystal turbulence 
and established that it is in the KPZ universality class. For flat initial conditions, their experimental asymptotic value 
for skewness was close to ~ 0.29 as suggested by their experimental plots. These values are displayed in Table 1 for 
comparison with our result. 


TABLE I: Values of Skewness in one dimension 


System of Study 

Reference 

Methodology 

Skewness 

SSM (flat) 

m 

Numerical 

0.28 + 0.04 

SSM (stationary) 

m 

Numerical 

Si 0.33 

DPRM (point-to-line) 

m 

Numerical 

-0.296 + 0.028 

DPRP (point-to-line) 

[26] 

Numerical 

-0.29 + 0.02 

RSOS (flat) 

[26,^ 

Numerical 

« -0.29 

TLC (flat) 


Experimental 

0.29 

PNG (curved) 

[20] 

Numerical 

0.2241 

PNG (flat) 

[20] 

Numerical 

0.2935 

PNG (stationary) 

[20] 

Numerical 

0.35941 

KPZ (mean fleld, flat) 

m 

Analytical 

±0.46 

Gombustion front (flat) 

[^ 

Experimental 

0.33 

Gombustion front (stationary) 

m 

Experimental 

0.32 

KPZ (present calculation) 

Eq. dZQl) 

Analytical 

0.3237 


We observe that our calculated result for the skewness is comparable with some of the experimental values and 
numerical simulations. It deviates from the non-stationary results and the deviation is more pronounced from the 
result for curved interfaces. This is expected as our calculations are applicable for a flat geometry in the stationary 
state. We also observe that the result of the mean field calculation deviates somewhat strongly from all other results. 
Since skewness is determined by the underlying probability distribution, its calculation following from the governing 
dynamics is of importance in inferring the universality class. Moreover, the existing studies indicate that the pdf is 
determined by not only the governing dynamics but also by the boundary conditions. Thus it may be said that there 
are different subclasses belonging to the same universality class. Different numerical values of skewness may thus be 
said to correspond to different universality subclasses although they may have the same scaling exponents for the 
correlation and response functions. 

Since the renormalization scheme involves calculations of the statistical moments in the large scale limit fc —>■ 0, such 
calculations are expected to lead to the statistical properties of the growth process at large scales. The fact that W 2 
and W 3 turn out to be infrared divergent implies a dominant role of the large scale fluctuations in determining these 
statistical moments. Moreover since the renormalization scheme involves calculations in the long time limit w —^ 0, 
such calculations are expected to capture the statistical properties of t he g rowth process in the stationary state. 

However, for a large system, achieving a stationary state is difficu lt j57l| . especially in experiments and numerical 
simulations, unless a stationary state is taken as an initial condition |58l |. To achieve a stationary pdf for the flat one 
dimensional KPZ problem, Imamura and Sasamoto took both sided Brownian motion as an initial condition [ 13,111 
and obtained the generating function for the replica partition function as a Fredholm determinant. This allowed for 
the calculation of the pdf which was found to approach the Baik-Rains Fq distribution in the long time limit. Krug 
et al. [ill investigated the stationary state skewness for the SSM model with random uncorrelated spins {ui = ±1). 
They obtained jS”! ~ 0.33 which agrees well with our calculated value. Maunuksela et al. [11 identified that the 
universality class of slow combustion fronts of a paper sheet belongs to the KPZ universality class on the basis of 
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scaling exponents. With the same experimental conditions, Miettinena et al. performed an experiment on paper 
burning to find the skewness from the pdf. They studied the height distribution of combustion fronts for flat initial 
conditions in the saturation regime and obtained the value S = 0.32 which agrees well with our calculated value. It 
may however be noted that Takeuchi [fi^ has suggested that their analysis of the pdf may need modifications. 

Although our calculated value for the skewness (0.3237) compares excellently well (within 1-2%) with the above 
mentioned stationary values (0.33 and 0.32), we observe that there is a slight departure (of about 10%) from the 
value S = 0.3594 coming from the PNG model [2^. This departure may be due to the dominating role of large scale 
fluctuations in determining the moments W 2 and W 3 . The infrared cutoff fj, may be interpreted as the inverse of 
the size L of the substrate and thus the calculations appear to be influenced by finite size effects. In spite of this 
slight departure, together with the agreements with the experimental results, our calculation for the skewness seems 
to identify the relevant universality subclass of the KPZ equation. 

Ideally speaking, full information about the pdf enables one to classify many seemingly similar problems into varying 
universality classes. However, an analytical calculation of the full pdf is an extremely difficult task, whereas the 
calculation of higher order moments such as the skewness is a more viable approach. Thus a classification scheme for 
universality beyond the scaling exponents can be formulated via the values of skewness for various processes. Takeuchi 
and Sano have proved through the TLC experiment that the KPZ class has geometry dependent subclasses in 
spite of having the same scaling exponents. Thus, to characterize the subclasses, knowledge of skewness and higher 
order moments is very much essential. The calculation of skewness, directly from the KPZ dynamics, is therefore an 
important step in identifying a universality subclass of the KPZ equation. 
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